Normal state thermodynamics of cuprate superconductors 
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We propose a microscopic explanation of the pseudogap features discovered in the normal state 
specific heat and magnetic susceptibility of cuprates. In the framework of the bipolaron theory of 
high-Tc superconductors we explain the magnitude of the carrier specific heat and susceptibility as 
well as their universal scaling with temperature over a wide range of doping. 



There is strong evidence for the normal state pseudogap in high-T c cuprates from magnetic susceptibility [Q, 
specific heat ||, angle-resolved photoemission (ARPES) j|, tunnelling and some kinetic measurements Q. One 
view supported by ARPES is that the gap reflects precursor superconducting correlations in the BCS-like state below 
some characteristic temperature T* without long range phase coherence ||. Testing of this hypothesis with specific 
heat [^| and tunnelling Q data, it is found that this view cannot be sustained. In particular, there is no sign that the 
gap closes at a given temperature T* , which rules out any role of superconducting phase or spin fluctuations Q . On 
the other hand, the strong-coupling extension of the BCS theory based on the multi-polaron perturbation technique 
firmly predicts the transition to a charged Bose liquid in the crossover region of the BCS coupling constant A ~ 1 [Q. 
The (bi)polaronic theory of carriers in cuprates, confirmed by infrared spectroscopy M and by the isotope effect on 



the carrier mass |£j], provides a natural microscopic explanation of the normal state gap |10|. Within the framework 
of the theory, the ground state of cuprates is a charged Bose-liquid of intersite bipolarons where single polarons exist 
only as excitations with an energy of A/2 or more. A characteristic temperature T* of the normal phase is a crossover 
temperature of the order of A/2 where the population of the upper polaronic band becomes comparable with the 
bipolaron density. Along this line the normal state kinetics of cuprates has been explained pd| , p^ | and a theory of 
tunnelling has been developed jl3| describing essential spectral features of STM and PCT conductance. In this Letter, 
we find a universal temperature scaling of the specific heat and magnetic susceptibility of Y Ba-iCu^O'j-i, and provide 
a microscopic explanation with bipolarons and thermally excited polarons. The central ideas of our model are as 
follows 

• Charge carriers are intersite real-space pairs of holes sitting on in-plane and apical oxygens |l4| ] . 

• In addition to introducing hole charge carriers, doping also introduces considerable disorder and localised states. 
Owing to interparticle Coulomb repulsion fHI , a localisation well contains either a bipolaron or an unpaired 
polaron but not both. The interaction of polarons and bipolarons in the extended states is taken into account 
within the Hartree-Fock approximation and included in their band dispersion. 

• At finite temperatures, a fraction of the carriers exist as unpaired hole polarons. These particles are responsible 
for the magnetic response of the system. 

We also employ the simplification that the tunnelling probability between localisation wells is negligible. This allows 
the partition function Zi for the localised part of the system to be written as 

zi = n Zi 

i 

Zf = l + 2e^ E ^ + e 2(M-Ei+A/2)/3 j (1) 

where we have assumed the no double occupancy condition. A, /i and E t are respectively the bipolaron binding 
energy, chemical potential and a single-particle energy level of the well, whilst [3 — 1/fcgT '. The point to note about 
equation ([j]) is that the localised partition function cannot be factorised into a product of two particle and one particle 
partition functions. The physics of localised bipolarons and polarons is thus not separable implying that only one 
density of states (DOS) profile should be taken for localised particles. The density of localised particles is determined 
by 
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with fli = —(3 1 log Z\. This gives 



ni=2 [ Pl {E)fi{E)dE (3) 
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where 

f l (E) = {l+g[0(E-p-A/2)}}- 1 (4) 

with g(£) — exp{£)cosh{^/2 + (3A/4)/cosh{^/2 — /3A/4). pi(E) refers to the density of localised states per spin. We 
can then write for the number conservation condition: 

2rib + 2n p + ri\ = x, (5) 

where nb, p is the density of delocalised bipolarons and polarons, respectively and x the doping per unit cell. For 
La2- x Sr x CuOi x is given by the atomic concentration of Sr whilst in Y Ba,2Cu?,0-!-&, x — 2(1 — <5)/3. The free 
particle density is given by 
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n b = / dEp b (E)f b (E) 
Jo 
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n p = / dEp p (E)f p (E) (6) 
Jo 

where fb(E) = {exp{(3(E — 2p — A)] — and f p {E) — {exp[P(E — p)] + l}" 1 , so that equation (5) allows us to 
determine the chemical potential p{T) if bipolaronic and polaronic DOS, pb yP (E) are known. The finite bipolaron 
bandwidth, the one-dimensional singularity of (bi)polaronic DOS |14|] , and a finite width of the localised tail give rise 
to a Shottky-like anomaly of the specific heat and a Curie- like temperature dependence of the susceptibility which 
are observed at high temperatures in overdoped samples as explained in Ref. |l2j . Here we consider the underdoped 
region, where the potential wells are deep and impurity-scattering broadening of the Van-Hove singularities (VHS) 
large due to ineffective screening by carriers. The previous analysis JlO| , p^[ showed that the characteristic width of the 
localised tails and VHS is above room temperature in underdoped samples. We can thus neglect any DOS structure 
for the relevant temperature range by taking P i(E) — Pp (E) — 2p b (E) ~ N(0) with iV(0) a single-particle DOS at the 
mobility edge, E = 0. The bipolaron chemical potential 2p + A is then pinned at the mobility edge, giving p — —A/2, 
as follows from Eq. (||) for fcsTiV(O) 1. This assumption greatly simplifies further calculations. Including the 
contribution of delocalised bipolarons, thermally excited polarons and localised carriers we obtain the total energy as 

E(T) = E Q + I" <%{Z[coth(0 - 1] + (4£ + /3A)[1 - tanh(£ + A/3/4)] + 2^[g(-^- 1 + l]" 1 }. (7) 
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Here Eq is a temperature independent (negative) constant. As a result we find a universal temperature scaling of 
the energy, E(T) = /(/3A), which allows us to extract the normal state gap from the experimental specific heat 
C = dE/dT without any fitting parameters as shown in Fig.l. In the low-temperature limit, f3A > 1 we get 
g(£) ~ exp(2£) and a linear specific heat with an exponential correction 



C ~ k B N(0)(3~ 
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■exp(-0A/2) 



(8) 



This result is in contrast with an expectation that the specific heat of nondegenerate bipolarons is temperature 
independent above T c . The random potential as well as a low-dimensional DOS pins the chemical potential at the 
mobility edge even in the normal state, so the bipolaron density (and hence the specific heat) is proportional to 
temperature. The latter leads to a temperature dependent Hall effect [ p"H and explains other anomalous kinetic 
properties of cuprates Jn|. Half of the bipolaron binding energy A/2, which is an energy gap between the bottoms of 
bipolaronic and polaronic bands has been estimated from 400^ to hOK depending on doping . In this temperature 
range one has to calculate E(T) and 7 = C/T by numerical integration of Eq.(7) with the result shown in Fig.l. 
There is a clear scaling of experimental 7 with /3A in a wide doping range of Y Ba2Cu^O-j-s- The corresponding 
values of A are shown in Fig. 2. They follow the same doping dependence as that determined phenomenologically 
and are of the same order of magnitude. It should be noted though that the d-wave approach taken by Loram et 
al gives consistently higher gap values than those found here. Nonetheless, d or s-wave like gaps in the DOS can be 
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obtained easily within this model by adjusting the nature of the particle-particle interaction or the dependence 
of the polaronic energy. Such discrepancies are not thus a significant problem. A drop of 7 at higher temperatures 
(Fig.l) is due to a finite bipolaron bandwidth as discussed above. 

Following Loram's analysis we compare 7 with the differential magnetic susceptibility x* — d{xT) / dT . The 
experimental data for x*{T) are perfectly consistent with our model. There are two contributions to the magnetic 
response, from delocalised (thermally excited) polarons, x p , and from localised ones, xi- For the first contribution we 
obtain by the use of the Kubo formula for free fermion magnetisation, 

Xp (T) = 2(j? B N(0)[exp((3A/2) + (9) 

where /is is the Bohr magneton. The single- well partition function in an external magnetic field, H is given by 

Zi = l + e 2(M--B+A/2)/3 + e (^-E+^ B H)P + e (p-E-^ B H)/3 qqs 

Differentiating twice the corresponding f2 potential with respect to the magnetic field yields 

Xi{T)=n%(3 [ dE Pl (E)f?(E), (11) 

where f[(E) = [1 + exp((3A/2)cosh{{E — \i — A/2)/3)] _1 is the distribution function of localised polarons. If DOS is 
a constant pi(E) = N(0), and temperature is low, j3A > 1 we obtain an exponential temperature dependence of the 
spin susceptibility as 

X (T) = x P (T) + X l(T) ~ 2/i|7V(0)(l + n/4)exp(-0A/2). (12) 

The numerical integration of Eq. (11) for the entire temperature range with the constant DOS yields a universal 
scaling of x* as a function of /3A. This is nicely confirmed by experiment, as shown in Fig. 3. It is remarkable, 
that with about the same A and DOS (see Fig. 2) one can describe both the specific heat and spin susceptibility of 
underdoped Y Ba2Cu^OT-s- This is at variance with some opinions that the experimental Wilson ratio is difficult 
to understand within the framework of our model. In fact, thermally excited polarons provide the spin susceptibility 
and a finite Wilson ratio close to the experimental one, while the binding energy of bipolarons is responsible for the 
normal state gap. We may therefore conclude that the formation of real space pairs (bipolarons) above T c and their 
partial localisation by the random potential are essential features in describing the normal state thermodynamics of 
Y BaiCuzOi-s and other cuprates exhibiting similar normal state gap. The bipolaron theory can explain such non- 
Fermi liquid features as a large carrier entropy, the gap above T c , temperature dependence of 7 and x an d their ratio. 
Another strong indication of the existence of bipolarons comes from the resistive and thermodynamic measurements 
in the critical region. A divergent upper critical field was measured in many cuprates as predicted by one of us [ jl8| , 
and the magnetic field dependence of the specific heat jump is just that of the charged Bose-gas p9[ . We greatly 
appreciate the enlightening discussions with A. Blackstead, J.R. Cooper, J.T. Devreese, J.D. Dow, A. Junod, H. 
Kamimura, W.Y. Liang, J. Loram, J.L. Tallon and G. Zhao. 
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Figure Captions 

Fig.l Universal scaling of j/k 2 B N(0) with 2k B T/A compared with theory (line) for YBa 2 Cu^0 7 -8 (N(0)= 1.17eV _1 
per spin). Fig. 2 Theoretical normal state gap as a function of doping. Fig. 3 Universal scaling of the differential spin 
susceptibility, x*(T)/n%N(0) = (xt xp - 0.39 x 10~ 4 emu/mole)//Lt3iV(0) compared with theory (line). 
For experimental data, see Loram et al pj\ 
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Figure 1: Universal scaling of j/kgN(0) with 2ksT/A compared with theory (line) for 

YBa 2 Cu 3 7 - 5 (N(0) = 1.17cV _1 per spin). 
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Figure 2: Theoretical normal state gap as a function of doping. 
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Figure 3: Universal scaling of the differential spin susceptibility, x*(T)/ 'fi^N(0) = (xt xp ~ 
0.39 x 10~ 4 cmu/mole)/^iV(0) compared with theory (line). 
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